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Abstract: The continuum version of the dual theory for a system of two-dimensional quan- 
tum disordered bosons with statistical particle-hole symmetry and at T = admits a new 
disordered critical point within the renormalization group calculation at fixed dimension. 
We obtain the universal conductivity and the critical exponents at the superfluid-Bose glass 
transition in the system: a c = 0.25(2e) 2 //i, v = 1.38 and z = 1.93, to the lowest-order in 
fixed-point values of the dual coupling constants. 



The problem of competition between superfluidity and localization in two dimensions 
(2D) has attracted lot of interest, as 2D is the lower critical dimension for both transitions 
when considered separately 0. Experiments on disordered thin superconducting films show 
a sharp separation between samples of varying thickness that exhibit diverging or vanishing 
conductivities as temperature is lowered, providing support for the concept of quantum 
(T = 0) superconductor-insulator (SI) critical point 0. Under the assumption that the 
amplitude of the superconducting order parameter remains finite at the critical value of the 
tuning parameter, it is sensible to describe the SI transition solely in terms of bosonic (phase) 
degrees of freedom in a random external potential. The resulting problem of disordered 
interacting bosons 0, Q is interesting on its own account, and has direct implications for 
the phase transitions in other systems like 4 He in porous media, Josephson junctions arrays 
or granular superconductors. Not a small part of the fascination with the SI transition in 
2D comes from a theoretical possibility that right at the transition point there is a finite dc 
conductivity, and moreover, that this conductivity in units of (2e) 2 /h should be a universal 
number [[[]. The universal metallic behavior at the SI transition is an option unique for 2D, 
which posed a challenge for analytical calculations of the critical conductivity , particularly 
in the disordered case where this quantity has hitherto been obtained only numerically ||. 
Generally, the critical behavior of disordered interacting bosons in more than one dimension 
has proven very difficult to address analytically ||, @, partly due to the fact that there 
is no easily identifiable correct mean-field starting point around which to set up a systematic 
perturbative calculation. 

In this Letter we present the results of a new approach to the quantum critical behavior 
of 2D disordered bosonic systems, applied to the paradigmatic problem of disordered lattice 
bosons at a commensurate density, interacting via short-range repulsion. Without disorder 
and at T = the system (described by the Hamiltonian (1)) undergoes a Mott insulator- 
superfluid (MI-SF) transition as the ratio t/U is varied, which is due to a commensurable 
boson density in the universality class of the 3D XY model J|]. With disorder the transition 
into a superfluid state is expected to always proceed from the gapless Bose glass (BG) phase 
@]. We study the SI transition in this system within the continuum version of the dual 
theory for the bosonic Hamiltonian fT0|| , which describes topological defects in the quantum 
ground state, and for the problem at hand takes the form of 3D classical anisotropic scalar 
electrodynamics (Ginzburg-Landau superconductor) in random magnetic field. The dual 
formulation of the problem may be thought of as a higher-dimensional analog of the Haldane's 
density representation of dirty bosons in ID [0. Since in dual representation the Gaussian 
point represents the interacting superfluid phase of the original bosons, one could expect 
that perturbing around it may be more sensible than in the original formulation, where 
one expands around the limit of non-interacting particles. Critical behavior of the isotropic 
version of the dual theory without disorder has been studied recently in different context 
by perturbative and non-perturbative methods JOJ], ||I2" |, ||13|| . While a small anisotropy at 



the critical point of the pure theory is irrelevant, the correlated random magnetic field is 
exactly marginal perturbation to the lowest order in disorder. However, we show that it 
generates a correlated random-mass term for the disorder field, which represents a strongly 



1 



relevant perturbation at the pure critical point. In one-loop approximation and for the 
physical initial values of the coupling constants the flow at the critical surface is attracted 
to a new disordered fixed point, which governs the critical behavior at the BG-SF transition 
in the system. We calculate the universal resistivity at the BG-SF transition, as well as 
the correlation length and the dynamical critical exponents v and z, as series in fixed-point 
values of the dual coupling constants, to the lowest non-trivial order. The comparison with 
previous studies and with experiments is discussed. 

Specifically, we are interested in the Hamiltonian for a collection of charge-2e bosons 



with short-range repulsion, U > 0, written in the number-phase representation jffl, |T0 



cos( 



where index i labels the sites of a 2D lattice, hi represents the deviations from a large 
integer average number of bosons per site, {hi} are Gaussian random variables with zero 
mean and width uih, and 0j is the phase variable canonically conjugate to the number of 
bosons, [4>i,hj] = i5ij. The Josephson term couples only the nearest neighbours. Using 
a set of exact and nearly-exact transformations on a lattice |TJj], |T5] Fisher and Lee \IU 



have demonstrated that the partition function for the disordered problem at T = may 
be written in the form of classical 3D (two spatial and one imaginary time dimensions) 
strongly anisotropic lattice superconductor in a random magnetic field. The dual lattice 
theory should be in the same universality class as the continuum field theory for a strongly 
type-II superconductor (A >> q 2 ) |T5[] : 

S = J d 2 fdz[\(V - iqA(r, z))*(r, z)\ 2 + /i 2 |^(f, z)\ 2 + ^(f, z)| 4 

+1(V x A(f,z)) 2 + |(V x A(f,z)) 2 + h(r)(V x A(r,z)) & ], (2) 

where the complex \l/ represents the disorder- field [0. The condensation of $ implies the 
proliferation of vortices in the ground state and consequently, the destruction of superfluidity 
by giving (Higgs) mass to the superfluids sound mode, represented by the dual gauge-field 
A (V • A = 0). Note that the random magnetic field along z (time) direction h(f) depends 
only on two (real space) coordinates, so that disorder is correlated along the third direction, 
as typical for quantum problems with quenched disorder. The condition on the average field 
h = corresponds to integer average number of particles per site. 

The SI transition in dual theory is reached by tuning the renormalized mass of the 
disorder-field to zero. Introducing replicas in standard way to average over disorder, the 
dual theory at m = may be parametrized by four dimensionless coupling constants: the 
dual charge q 2 = q 2 /p, the quartic term coupling A = X/p, the anisotropy parameter a, 
and the coupling Wh = Wh/p, which appears in the disorder-induced term of the form 
— (wh/2) Y^a8 /(V x A a (r, z))z(V x Ap(r,z'))$. p is an arbitrary infrared scale. Besides 
these, the random magnetic field generates yet another quartic coupling constant w = wjisp 2 
via the term —(w/2) Y, a a I l^a(^*; ^OPl^M^ z ')\ 2 i n the replicated theory, with a meaning of 



2 



width of the Gaussian distribution of the random-mass term (again correlated along z axis) 
for the disorder-field. To determine the flow of the five coupling constants we performed a 
perturbative RG calculation in fixed dimension ||T3] ; [13|, since we will eventually be 



interested in the critical conductivity, which is finite only in 2D. The gist of the method is to 
express the universal quantities as series in renormalized, instead of bare, coupling constants, 
so they tend to finite values in the infrared limit, determined only by the location of the 
attractive fixed point. We define the renormalized couplings A and w as the corresponding 
quartic vertices right at m = and at the usual symmetric point in the momentum space: 
ki ■ kj = (4:5i 7 j — l)p 2 /4, i,j = 1, 2, 3, with additional condition ki jZ = k 3tZ = —ki tX = —k^ z to 
accommodate the correlated nature of the disorder vertex w along z-axis. The polarization 
diagram which renormalizes q, a and Wh is evaluated at the external momentum c • p, where 
the parameter c will be specified shortly. The /3-functions for the anisotropy a and the 
random-field width are: 

^ = -«(^ 2 + 0(g 4 )) + 0(a 2 ), (3) 

^ = MI-Va) + 0(wI), (4) 

where t = ln(l/p) — > oo, in the infrared limit, and renormalized couplings are assumed. 
Eq. 3 is a one- loop result, and suggests that a small anisotropy is irrelevant at any fixed 
point with a finite dual charge. Although the dual theory should be strongly and not weakly 
anisotropic [IIJ, we do not expect stable fixed points with a ^ 0. This is known to be 
true in the classical scalar electrodynamics close to four dimensions and for a large number 
of complex field components JT7|. In contrast to eq. 3, eq. 4 is exact to the lowest order 
in Wh- This follows from the fact that the transverse part of polarization is diagonal in 
replica indices to the lowest order in Wh, so that the renormalization of Wh to this order 
entirely comes from rescaling of the gauge-field. The anomalous dimension of the gauge-field 
propagator is t/a = 1 in 3D at any fixed point with a finite dual charge |12| . In particular, the 
linear term in eq. 4 vanishes at the pure MI-SF fixed point (see below). This observation may 
explain why numerical calculations at commensurate boson densities always have difficulties 
in detecting the intervening Bose glass phase at weak disorder |18 |. It also seems consistent 
with the result of Singh and Rokshar || , where a weak disorder in bosonic Hubbard model 
is found to be irrelevant, but only unusually weakly so. Fate of the coupling constant Wh 
is determined by the higher order terms in eq. 4. To proceed, we will assume that Wh is 
irrelevant at the MI-SF fixed point of the theory fl9 |, and that its only role was to generate 
a strongly relevant coupling w, before renormalizing to zero. Hereafter we will therefore set 
both a = Wh = 0. 

One-loop /3-functions for the remaining three coupling constants are then: 

^ = g -16 9 ' (5) 
f = A(l + \q 2 + 2(V21n(l + V2) + ± ln(V2 + Vs))w) - ^|±1a 2 - ^=q\ (6) 
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f = tf,(2 + \? - ^ + 2( 71 ln(1 + ^ + 7s la{V ~ 2 + ^ (7) 

For w — these equations reduce to those studied previously in the context of critical 
behavior of superconductors fl2fl . In that case for a choice c > 5.17 |20| there are four fixed 
points in the theory: Gaussian (g 2 = A = 0) and 3D XY (g 2 = 0, A xy = 2.09), both unstable 
in direction of the dual charge, tricritical (g 2 = 16/c, A = A_), unstable in A-direction, and 
the MI-SF critical point (q 2 = 16/c, A = A+), believed to be of "inverted" XY type @ (see 
Figure 2. in ref. 12). The Gaussian and the tricritical fixed points are connected by the 
straight separatrix determined by the tricritical value of the Ginzburg-Landau parameter 



k = y\/2q 2 = k c in the problem. On the other hand, the tricritical value of k has been 
estimated both analytically (TB| and numerically pi|, which allows us to fix the parameter 



c to match that number: n c = 0.42/a/2 requires c = 20. It is worth noting that for w — 
our calculation gives reasonable estimates of the critical exponents both at the unstable XY 
(u = 0.63) and at the stable MI-SF fixed point (u = 0.61, rj = -0.20). 

A small w is relevant at the MI-SF fixed point, as follows from eq. 7 and from the Harris 
criterion [22|. Both perturbative and non-perturbative |ll|, [II], [IS] calculations of the 



correlation length exponent at the pure critical point yield values v < 1, suggesting relevance 
of the correlated random-mass term. Starting from the type-II region (A >> q 2 ), the flows 
are attracted by the disordered critical point at q 2 = 0.8, w c = 3.71 and A c = 29.72, which 
we identify as the BG-SF critical point. The BG-SF critical point exists for any choice of 
renormalization procedure (parameterized by c), unlike the charged fixed points in the pure 
theory which exist only for c > 5.17. This is akin to the classical scalar electrodynamics 
close to 4D [p3"| , where randomness restores the critical point in the theory. To the lowest 
order in the critical coupling constants the exponents at the BG-SF fixed point are: 

, = l(l + ^-^- 4 ^ ) = 1.38, (8) 

z = i + -^ = im. (9) 

4 v ' 

Note that v > 1, as expected |p2 |. The dynamical exponent is very close to two, which was 



conjectured to be an exact result M in 2D. The exponents are also in good agreement with 
Monte Carlo calculations of Wallin et al. || [y = 0.9 ±0.1, z = 2 ±0.1) and real-space study 
of Zhang and Ma [| (u = 1.4, z = 1.7). 

The disorder-field \l/ describes condensation of topological defects in the bosonic ground 
state, and duality of charges and vortices implies that the dimensionless conductance for 
vortices is the resistance for the original bosons. Kubo formula for the conductivity (or 
conductance, in 2D) of the vortices described by the dual theory then leads to the critical 
dc resistivity at T = 0: 

R c = ^% lim±[2<¥*(0)¥(0)) - / d 2 rdze^(j x (r,z)j x (0,0))], (10) 
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where j is the current of the dual field. To the lowest order we obtain 

,7T 7T ~ 2.94 „ . h 
8 16 7r (Se)" 1 

The first term describes the Gaussian part of the resistivity at the critical point f|, while 
the second term accounts for the lowest order wave-function renormalization. The third, 
disorder term, is the main source of dissipation. Inserting the values of the critical coupling 
constants into eq. 11 we estimate the universal conductivity for the bosons: 

o c = R- 1 = 0.25^. (12) 

The result is somewhat larger than in Monte Carlo calculations of Wallin et al. ||: cr c = 
[0.14 ± 0.03](2e) 2 //i, and smaller than in Batrouni et al. @: a c = [0.45 ± 0.07](2e) 2 //i. Our 
preliminary calculations suggest that the higher order corrections tend to lower the value of 
a c . Remark however that the value of a c is already smaller than in the pure case ||. 

The reader has certainly noted that the BG-SF fixed point in question is a strong- 
coupling one, and the numerical values of critical quantities cited above need to be taken 
with some reservation. There is no small parameter in the problem and we had to rely on 
a less controlled calculational scheme in fixed dimension [|TT|. Also, a different tricritical 
value of GL parameter k c in zero-disorder scalar electrodynamics would dictate a different 
choice of the parameter c, and would affect the position of the SF-BG fixed point within 
our calculation. The dependence of the critical quantities on this variation seems weak, 
fortunately. For instance, assuming almost a twice larger value for k c = 0.8/ \/2 [[ll| we 
would find v = 1.55, z = 2.14 and a c = 0.18(2e) 2 //i. Our calculation allows for a systematic 
improvement by going to higher orders in the dual coupling constants and using standard 
resummation techniques [pl| . Although the perturbation series in eqs. 3-7 are likely to be 
only asymptotic (but Borel summable [|4|), it is encouraging that the lowest order calculation 
already produces sensible results. 

We have observed that our numbers for the critical quantities are not too different from 
those obtained by other methods performed on the Hamiltonian (1) || or on its variant 
||, but at a non- commensurate density of bosons. Without disorder, at non-commensurate 
filling factors the MI-SF transition is in a different universality class, and it is in fact mean- 
field in character ||. One would expect on physical grounds that commensurability may be 
irrelevant at the disordered critical point || . Deviation from a commensurate density would 
correspond to a finite average magnetic field h in dual theory, whose scaling dimension at 
the BG-SF fixed point we do not presently know. We have argued however, that even at 
commensurate densities the transition is always controlled by the disordered BG-SF fixed 
point, due to the generated random- mass term for the disorder field. 

Comparison with experiments is made difficult by the fact that different measurements 
lead to rather different values of the three critical quantities we considered 0. We may 
still note that our numbers for the critical exponents fall within the range vz = 2.8 ± 0.4 
obtained from the scaling of resistivity on the insulating side of the SI transition in thin Bi 
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films |2[. The value of the critical conductivity in thin films differs by a factor of ~ 4 in 
various experiments |§, and its universality may be questioned. Many experiments show 
er c ~ (2e) 2 /h [fj], somewhat larger than our estimate. Note also that experimentally one 
typically measures u = 0, T — > conductivity f7j, which is the opposite order of limits than 



assumed here. It has been argued recently that the two limits may not commute [25], and 
that the critical conductivity has a non-trivial structure as a function of u/T. 

In conclusion, we considered the critical behavior at the superfluid-insulator transition in 
a system of 2D disordered bosons with commensurate density and at T = 0, by studying the 
continuum version of the dual theory that describes vortices in the superfluid ground state. 
Although disorder, which initially appears in the dual theory as random magnetic field, is 
argued to be marginally irrelevant at the critical point of the pure theory, it generates another 
strongly relevant disorder-like coupling constant, and the transition is ultimately governed by 
the BG-SF critical point. The universal resistivity at the BG-SF transition and the critical 
exponents have been expressed as perturbation series in the critical coupling constants, and 
the lowest order results compare well with the Monte Carlo calculations performed on the 
original boson Hamiltonian. 
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discussions. This work has been supported by NSERC of Canada and Izaak Walton Killam 
foundation. 
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